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Abstract. We derive black hole masses for a sample of about 300 AGNs in the redshift range < z < 2.5. We use 
the same virial velocity measure (FWHM H/3bc) for all sources which represents a significant improvement over 
previous studies. We review methods and caveats for determining AGN black hole masses via the virial assumption 
for motions in the gas producing low ionization broad emission lines. We derive a corrected FWHM(H/3bc) measure 
for the broad component of H/3 that better estimates the virialized line emitting component by comparing our 
FWHM measures with a sample of reverberated sources with H/3 radial velocity dispersion measures. We also 
consider the FWHM of the FellA4570 blend as a potential alternative velocity estimator. We find a range of black 
hole mass between logMsH ~ 6.0-10.0, where Mbh is in solar masses. Estimates using corrected FWHM(H/3), 
as well as FWHM(Fell) measures, reduce the number of sources with log Mbh > 9.5 and suggest that extremely 
large Mbh values (log A/bhJJ 10) may not be realistic. Derived L/LecM values show no evidence for a significant 
population of super-Eddington radiators especially after correction is made for sources with extreme orientation 
to our line of sight. Sources with FWHM(H/3bc) ^ 4000 km s~^ show systematically higher L/LEdd and lower 
Mbh values than broader lined AGNs (including almost all radio- loud sources). 

Key words, quasars: emission lines - quasars: general - line: profiles - black hole physics 



1. Introduction 

Gravitational accretion onto supermassive black holes is 
generally accepted as the ultimate energy source of Active 
Galactic Nuclei (AGNs). The last decade has seen a major 
effort to derive reasonable estimates of black hole masses 
(A^bh) by assuming virialized motions in the broad line 
emitting gas: 



A^BH 



G 



(1) 



where v is the velocity dispersion of the emit- 
ting gas at distance tblr- The velocity dispersion 
can be written as u = /FWHM where FWHM is 
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the full width half maximum measured for a suit- 
able emission line. The factor / depen ds on the ge- 
ome try and detai l s of the kinem atics f Krolikl l20nit 
McLure fc Dunlonl l200lt lOnken et a l.. 2004. The use- 
fulness of the virial assumption for Mbh determina- 
tion is best exploited using reverberation-mapping s tud- 
ies, e s pecially for the Balmer line s (Koratka r fc GaskellL 
| l99lt IPeterson & Wandel Il999t Iwandcl et all Il99a 
iKaspj et al.'. l 200(i IPeterson et aU i2004^) . iKaspi et all 
(|200(i, .200511 derive a relation between Broad Line Region 
(BLR) distance tblr and continuum luminosity 



''BLR oc {XLxY 



(2) 



where the exponent a is constrained between 0.5 and 1.0, 
and is most likely k, 0.6-0.7 for broad H/3 (H/3bc) and 
optical continuum luminosity. 
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There are caveats associated with this method. The 
virial assumption is not hkely to be generally valid for the 
emission line gas in AGNs. It has been known for sev- 
eral decades that different emission lines in a source can 
show different width and profile shape fe.g.. lde Robertisl 
ll985HSulenticlll989t) . The virial assumption implies that 
the velocity dispersion will steadily decrease with dis- 
tance from the central black hole oc ' . The time lag 
between continuum fluctuations and corresponding emis- 
sion line responses will therefore anti-correlate with line 
width. This trend has bee n confirmed in a few objects 
l|Peterson fc Wandej l200o|) . 

Observations however show that profile width and 
shape depend on the ionization potential. The strongest 
high-ionization lines (HILs; e.g. CivA1549) often display 
blueward asymmetric profiles, or even centroid blueshifts 
(up to several 1000 km s~^) with respect t o the best esti 
mates of the rest frame of the source (e.g 
MaTzia,ni et aHll 99fiHRirha,rds et a1.Ll2002t 



Gaskell'. 11982 ; 
Bachev et al 



2004) . Bhieshifts like those observed for CivA1549 are in- 



dicative of obscuration and radial motions which invali- 
date the virial assumption. 

Low ionization lines (LILs) like H/3 are the best can- 
didates for emission arising from a virialized medium. All 
or much of the Balmer (and Fell) line emission is thought 
to arise from an accretion di sk or a flattened cloud dis- 
tribution near the di sk (e.g., ICollin-Souffrin et all Il988t 
iMarziani et alJ . fl996(l . A major caveat, especially for the 
Balmer lines, involves the possibility that: (1) there are 
two or more emission components in a line, and (2) only 
one of them may arise in a region where the virial assump- 
tion holds. The Balmer lines do not usually show very 
large line shifts (i.e., shift < FWHM) although profiles can 
be very asymmetric. Both red and blue asymmetries are 
observed for H/3bc which is the most studied line because 
it is relatively unblended a nd is observable with optical 
spectrometers up to z » 1 llOsterbrock fc Shudeii Il982 ; 
SulenticLll989HStirnel.ll99(llSulentic et allll99flHCorbin . 



199lUCorbin fc Francij. Il994l) . 



The most ambiguous sources from the point of view of 
Mbh determination sho w FWHM(H/3Rn ) ^ 4000 km s^^ 



(Population B, following Sul entic et al.L l2000bt) and red- 
shifted profiles and/or red asymmetries. Not all parts of 
the Hz? profile respond to continuum changes in the same 
way implying that some of the line emitting gas may be 
optically thin or, less likely, is not exposed to the variable 
Hi ionizing radiation (we will return to this issue in H5.2(l . 
The broad line profile in Pop. B sources may be due to two 
distinct emitting regions: (1) an optically thick classical 
BLR and (2) a broader and redshifted very broad com- 
ponent that may be optically thin or marginally optically 
thick to the Lyman continuum, ori ginating in a distinct 
Very Broad Line Region (VBLR.) jMa rzi^ni fc Sulentic] . 
ll99.lHShield~ et al. l ll99fiHSulentic et al. U2000(j) . The red- 
shift of the VBLR component raises doubts that it arises 
from virialized gas. A strong BLR response to continuum 
fiuctuations, coupled with a weak or absent response of the 
VBLR component, can lead to an overestimate of FWHM 



for the virialized BLR compon e nt resulting in an o veres- 
timate of MpTT llWandel et all Il999t iKaspi et all I2OO0I: 
IVestergaardl 120021 

Sources with FWHM(H/3Rn) ^ 4000 km s"! 
(Population A. ISulentic et allboOOhl) should provide more 
reliable Mbh estimates. The H /^bq profile is usual l y well 
fit with a syrnmetric function <)Veron-Cettv et al~l l200lt 
ISulentic et all l2002f) and the BLR emission is thought 
to arise from a Keplerian disk. The most unreliable 
Pop. A sources, in a disk emission scenario, should be 
those observed near face-on where the rotational (i.e. 
virial) contribution to FWHM(H/3bc) is minimal. At least 
some of the face-on sources may be identified as the so- 
called "blue outliers" which show a weak and significanth 
blueshifted fOi i il AA4959,5007 fin es ("Zamanov et al.', 2C 
iMarziani et^aP. l2008bt lAoki et al.„ .2005; .BorosonL 
The inferences in this, and the preceding paragraph 
have emerged from the Eigenvector 1 (El ) scenario that 
we have pursued for the past 5-|- years ISulentic et all 
l2000albL following iBoroson fc l(l992f) ). The results 

reported in this paper appear to confirm them. 

There are additional caveats connected with us- 
ing FWHM(H/3bc) for Mbh derivations. Line asym- 
metries can affect Mbh estimates for both Pop. A 
and B sources. The Vr displacement of the line cen- 
troid at fractional intensities (at half maximum, c(^)) 
can give useful information about the uncertainty of 
FWHM(H/3bc) measurements. The average displacement 
value is a fe w hundred kin s^J_with median c{\) « 
400 km s~i l|Marziani et all l2003a|) . This implies that 
deviations from a symmetric profile can affect Mbh 
estimates by « 40 %. Uncertainties are also intro- 
duced by: (1) contamination from overlapping/nearby 
lines such as Fell, H eiiA4686 and fQii i1AA4959 5007 
llOsterbrock fc Sbudeil 11981 Ide R.obertij. EM M 
Il988l: l.lackson et all Il99lh . (2) FWHM measures based 
on single-epoch observations, (3) low S/N spectra and (4) 
spectra without H/3 narrow component (H/3nc) subtrac- 
tion. The tblr - L\ relation is also not free from un- 
certainties. Reverberation mapping-based A/bh determi- 
nations are certainly affected by the non-negligible radial 
extent of the optically thick BLR. The derived tblr is 
not a very well defined quantity and a is also somewhat 
uncertain because of the intrinsic scatter in the correla- 
tion. Finally, reverberation data does not exist for high 
luminosity/rcdshift quasars requiring an extrapolation of 
the tblr - L\ relation in order to estimate Mbh for these 
sources. 

Uncertainties for A/bh derivations using single profile 
observation of H/3bc are estimated to be a factor of 2-3 
(at a 1 fj confidence level), but may be as low as 30% if 
the velocity dispersion of the varia ble part of H/^bc profil e 
is employed as a virial estimator ijPeterson et alll2004D . 
Mbh estimates based on the virial relation retain a statis- 
tical validity considering that AGNs span a 5dex range in 
A/bh (10^ Mo< A/bh ^ 10^" Mq) and that the relation 
has now been applied to large samples of objects (^ 10'^; 
iMcLure fc Duulon, 2OO4I). 
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With these considerations in mind, and supported by 
previous results, we use the virial relation to compute 
Mbh for « 280 AGNs in our El sample of low red- 
shift (0 ^ z ^ 1) sources (jMarziani et all l2003b|) sup- 
plemented with 25 intermediate redshift/high luminosity 
{1 ^ z ^ 2.5) quasars. We derive Mbh and L/LEdd values 
in several ways. New VLT-ISAAC data are presented for 9 
sources whi ch supplement data a lready published for 
sixteen quasars (|Sulentic et alll20fl4l) . Line measures are 
presented in We use H/3bc and FeiiA4570 line widths 
in a consistent way over the redshift range ^ z ^ 2.5 (a 
range of 10^ in luminosity). We compute black hole mass 
Mbh, Eddington ratio L/Lsdd (iSJ and we discuss how 
mass determinations might be improved f ii5.2|l so that the 
evolution of Mbh and L/ L^^^ with redshift can be con- 
sidered (SO and 



2. Observations and Data Reduction 

New intermediate redshift data were obtained between 
10/2003-03/2004 in service mode with the infrared spec- 
trometer ISAAC mounted on VLTl (ANTU) at the 
European Southern Observatory. Each spectrum corre- 
sponds to a wavelength range (IR windows sZ, J, sH) that 
covers the region of redshifted H/3bc and FeiiA4570 or 
FeiiA5130 at least in part. Reduction of quasar spectra 
and standar d stars followed exac tly the same procedures 
described in lSulentic et al.1 l(2004l) . Wavelength calibration 
yielded rms residuals of 0.4, 0.6 and 0.9 A in the sZ, J and 
sH windows, respectively. Absolute flux scales of the spec- 
tra will be inaccurate because atmospheric seeing almost 
always exceeded the slit width (w 0".6) resulting in sig- 
nificant light loss. 

Table n summarizes the new observations and the ba- 
sic format is given below the table. All sources come from 
the Hamburg-ESO (HE) quasar survey, which is a flux 
limited (with lirniting ttir « 17.5), color-selected survey 
i Wisotzki erall 12000 h . Column 2 of Table [T] hsts the 
blue apparent magnitud es from the HE s urvey papers 
l)Reimers et all [1996; Wisotzki et all 12000 1) while Col. 3 
lists the source redshift z computed as described in section 
121 The brightest sources of the HE at intermediate red- 
shift were preferentially selected. Col. 4 indicates whether 
[Olll]A5007 was used to compute z as indicative of the 
source rest frame. The absolute magnitude Afs reported in 
Col. 5 was computed by assuming Hq — 70 km Mpc"^ , 
and relative energy density ilM = 0.3 and JIa = 0.7. The 
k correction was computed for a spectral index a = 0.6 
{S oc j/~°). Col. 6 gives the ratio of log specific fluxes at 6 
cm and 4400 A (logi?k). In most cases only NRAO VLA 
Sky Survey (NVSS) upper limits are available. Columns 
7-12 give details of the observations explained in the foot- 
notes. The continuum S/N estimate given in Col. 12 was 
measured using a small region of the spectrum that was 
as flat and free of line emission as possible. 



3. Data Analysis 

3.1. Redshift Determination and Rest Frame 
Corrections 

Small offsets are present in the wavelength calibration, 
because the arc lamp frames were obtained in daytime, 
and therefore usually after grism movement. A correc- 
tion for these shifts was obtained by measuring the cen- 
troids of 2-3 OH sky lines against the arc calibration 
and calculating the average difference, which reached at 
most 6.5 A or 2.5 pixels in either direction. Rest frame 
determination for the 9 new sources was usually esti- 
mated from the H/3 peak redshift (assumed rest frame A — 
4861.33 A). [Oiii]A5007 (assumed rest frame A = 5006.85 
A) yielded a consistent measurement in only two sources 
(where the results were averaged. Two sources show no 
clear detection of [Oiii] AA4959,5007 while the remaining 
five show a significant disagreement between [Oiii]A5007 
and H/3nc- In these sources (we call the extreme exam- 
ples of them "blue outhers") it is not advisable to use 
[Qui ] A5007 for redshift det e rmination dZamanov et al' , 
2002t iMarziani et all l2003bt lAoki et al. I l2005t iBoro'sml 



200,51) . The adopted estimate was used to deredshift the 



cor- 



spectra while the dopcor IRAF task applied a (l + z)^ 
rection to convert observed specific fluxes into rest frame 
values. Fig. n shows the deredshifted spectra. 



3.2. Continuum and Feu subtraction 

Our spectral analysis made use of standard IRAF tasks 
with the first step involving continuum modelling and sub- 
traction. Using Chebyshcv polynomials of 

3rd 2nd 

der, a reasonably smooth continuum subtraction was ob- 
tained for all sources. To estimate errors in the continuum 
assessment introduced by noise, we also defined a mini- 
mum and a maximum continuum. Continuum fluxes were 
chosen at about -3(t (minimum) and +3a (maximum), 
where a is the noise standard deviation from the most 
likely continuum choice. Errors in continuum placement 
defined by difference between the extreme continua and 
the most probable one were then propagated according 
to standard error theory. The results of this procedure 
are cons istent with continuum fits employing very si mple 
models (" Malkan fc SargentL Il98l 'Shang et alll2004j) . We 
assumed that the continuum underlying the H/J spectral 
regions is due to two components: either a blackbody of 
temperature 25000° K or a power law of slope b = 0.7 
(/^ cx assumed to be valid only locally around H/3). 
In 5 sources the sole black body component produces a 
good fit; in 2 the blackbody component is dominant, and 
only in the remaining 2 cases the power-law alone can pro- 
vide a good fit. This method has some limitations due to 
the small spectral bandwidth covered by our spectra, to 
the relative strength of Fell and to internal reddening ef- 
fects. Since we did not attempt to change the blackbody 
temperature nor the slope of the power- law, we adopt the 
empirical continuum which is visually more accurate. 
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The emission blends of Fell were subtracted using 
th e template method based upon the spec t rum of I Zw 
1 l|Boroson fc GreenL Il992t iMarziani etall l2003a|) . The 
strongest Fell blends fall in the wavelenghth ranges 4450- 
4600 A (blue blend: FellA4570) and 5200-5600 A (red 
blend: FeiiA5130). The method includes the following 
steps: (1) the template intensity is scaled to roughly match 
the observed spectrum, (2) a Gaussian broadening factor 
is estimated from FWHM H/3bc (3) a careful estimate 
of minimum and maximum plausible broadening factor 
is made to set a reliable ±3ct uncertainty and (4) the 
template intensity is adjusted as necessary after broad- 
ening. The resultant template was subtracted from the 
continuum-subtracted spectra. The blue side of the spec- 
trum including FeiiA4570 is missing, or not fully covered, 
in several sources (HE 0946-0500, HE 1003-1-0149, HE 
1017-0009, HE 1249-0648 and HE 1258-0823). In these 
cases the best template fit was achieved for the red blend 
and a fixed ratio between the red and blue blends was as- 
sumed to reproduce FeilA4570. Fig.Elshows the estimated 
Fell emission (green lines). Fell was subtracted before con- 
tinuum fitting in sources with strong Fell emission. 

The [Oiii]AA4959,5007 lines were measured after Fell 
subtraction and taking into account the following con- 
siderations: (1) the flux ratio between [Oiii]A5007 and 
[Oiii]A4959 should be «3, (2) both lines should show 
identical profiles and (3) any H/3bc emission underlying 
[Oiii] AA4959,5007 is expected to have a smooth shape. 
Sources HE 1249-0648 and HE 1258-0823 see (Fig. 0) 
show a small bump at A « 5016 A which is either an Fell 
subtraction residual or emission from HeiA5016 (weak red- 
shifted [Oiii] A5007 would be almo st unprecedented). Afte r 
subtraction of H/3nc (following IMarziani et al.L l2003a|) , 
the H/3bc profile was fit with a high order spline function 
(IRAF task SFIT). This procedure does not yield a model 
fit but only an empirical fit that smooths the noise and 
reproduces the main features and inflections in H/3bc- 



4. Immediate Results 

4.1. Line Measurements and Uncertainties 

Line Fluxes and Equivalent Widths Table|31gives line mea- 
surements for the new VLT spectra with the basic for- 
mat given in the footnote. Cols. 3 and 4 give equivalent 
width (EW) measures for H/3bc and FeiiA4570 respec- 
tively. We evaluated uncertainties associated with the con- 
tinuum level (derived from the minimum and maximum 
reasonable continuum estimates) and line flux errors es- 
timated from the S/N . These estimates were combined 
quadratically to obtain uncertainties for EW measures. A 
similar procedure was applied to obtain uncertainty esti- 
mates for H/3nc and [Oiii] AA4959,5007. The EW uncer- 
tainty for H/3nc was derived from the estimated maxi- 
mum and minimum possible H/3nc component in the H/3 
line. The relative error of the H/3nc Aux can be large (see 
Table 121) and in some cases an H/3nc component may not 
be present. 



Fe iiopt width Fell emission is heavily blended so that 
widths of individual lines must be obtained from the 
best broadening parameter that was used for the tem- 
plate fit. This requires that we assume a constant width 
for all Fell lines which so far appears to be reason- 
able. FWHM(FeiiA4570) values derived from the tem- 
plate broadening factor are reported in Col. 5 of Table 
13 Uncertainty estimates for FWHM(FeliA4570) were ob- 
tained by increasing/decreasing the broadening factor un- 
til we could detect significant changes in the best fit. 
Simulated data reveal that it is possible to estimate the 
FeiiA4570 width up to FWHM(FeiiA4570) « 6000 km s^^. 
Due to the very large uncertainty of FeiiA4570 width de- 
termination for individual sources, Mbh estimates based 
on FWHM(FeiiA4570) measures are used mainly for con- 
firmatory purposes of statistical trends detected with H/3. 



H/3bc Line Profiles Measurements of FWHM(H/3) to- 
gether with other important line parameters like asymme- 
try index, kurtosis and line centroid at various fractional 
intensities were derived using a FORTRAN program devel- 
oped for that purpose. These parai neters are the same as 
defined in sev eral pre vious papers l|Marziani et al . ', '1996", 
|2003a; Sulent ic et all E'0041 and are given in TableEl Each 
line measure is followed in the next Col. by its appropri- 
ate uncertainty. Cols. 2 and 4 give the Full Width at Zero 
Intensity (FWZI) and FWH M. Col. 6 give s the a symme- 
try index (Al) as defined in lSulentic et all l|2004l) . Col. 8 
lists kurtosis values while the remaining part of Table 01 
lists measures of the H/3bc centroid at various fractional 
intensities, (in km s^^). All uncertainties represent the 2cr 
confidence level. 

The dichotomy in H/3bc profile shape ( and many other 
prope rties) between Population A and B l|Sulentic et all 
l2002l) is seen in the new source measures and in th e 
rest of our higher redshift sample JSulentic et all 120041) . 
Redward asymmetries (A.I. > 0.2) are most often found 
in Pop. B sources. They also show H/3bc profiles that 
are best fit with Gaussian functions, and some profile ap- 
pear composite. A few sources appear to deviate from the 
trend found in previous work t hat FWHM of H/^Rr: an d 
the Fell lines are very similar l|Marziani et al.l l2003al(J) . 
HE 1249-0648 and HE 1258-0823 show FWHM(H/3bc) 
> FWHM(FeiiA4570). While FWHM(FeiiA4570) is sub- 
ject to large uncertainty, the difference is confirmed by 
careful reinspection of these spe ctra. This condition is 
seen in only 2/215 sources in the IMarziani et ^ 1 2003ah 
sample. Both (IRAS 07598-1-6508 and Mki i 235) are BAL 
QSO s which are also FIR ultra-luminous ijSulentic et all 
l2006i) . In addition to a CivA1549 BAL with high ter- 
minal velocity, these objects have a strong, blueshifted 
CivA1549 emission line component. Blueshifted Balmer 
emission is probably associated with the high-ionization 
gas emitting CivA1549. The absence of any detectable 
[Oiii] AA4959, 5007 emission, along with a possible Hei fea- 
ture at A5016 A support the possibility that HE 1249-0648 
and HE 1258-0823 could be BAL QSOs similar to IRAS 
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07598+6508 and Mkn 235. Further support comes from 
an inspection of the HE survey spectra of both objects 
which indeed show blueshifted broad absorption in the 
MgiiA2800 hne. 

5. Black Hole Mass and Eddington Ratio 

5.1. Basic Equations 

One can write the velocity dispersion v in Eq. ^ as w ~ 
73/2 FWHM(H/3bc) in the case of randomly oriented ve- 
locities projected along the line-of-sight. The expression 
for black hole mass is then: 



A^BH 



3rBLRFWHM(H/3Bc) 
4 



G 



(3) 



In the absence of reverberation data one must rely on 
the correlation between reverberation radius and source 
luminosity tblr c< (L.=^inn )"i where L^^nn is the s pecific 
luminosity at A « 5100 A |Kasni et airi200nL l2005l) . This 



formula relates BLR distance from the central continuum 
source an d specific lurninosit y (ergs s^^A"^) near 5100 A. 
Following iKaspi et alJ l)2005|) one can write: 



r^LR. ~ 0.697 • 10 



17 



ALa (5100A) 
10^4 ergs-i 



0.67 



(4) 



We can make two estimates of A/bh using these rela- 
tionships assuming two different values for XL\. The first 
method derives A^b from tabulated values of y or i? ap- 
parent magnitudes: 



ALa (5100^ 



3.14 • io35-0-4(Mb) 



ergs s 



(5) 



A second possibility involves estimating the specific 
luminosity near 5100 A directly from flux measures using 
the adopted continuum fits for our spectra. More explicitly 
XLx is computed as follows: 



XLx (5100i 



47r4A/A(5100A) 



ergs s 



(6) 



where dp is the redshift derived distance and fx is the 
specific flux in the rest frame at 5100 A (after correction 
for Galactic extinction ^b at the observed wavelength). 
The second choice has the advantage that the continuum 
fx and FWIIM(H/3bc) are measured from the same spec- 
trum. Expected light losses are estimated to be « 35 % 

of the quasar flux with average Paranal seeing (0."6). We (c) profile asymmetries frequently observed in Pop. 

apply this average correction to the observed flux before 
computing A/bh- 

Substituting the expression for tblr, in the mass for- 
mula one obtains the following relation: 



reported in Table |21 AH sources with z ^ 0.8 are from 
the VLT sampl e. We do not show Af RH values derived 
from M-Q as in iMarziani et al.l ( 2003lJ l for clarity, since 
they basically confirm the same trends obtained from the 
specific fluxes. 



5.2. Improved Mbh Estimators? 

As briefly summarized in there is some evidence sug- 
gesting that the virial assumption is reasonable for LILs 
in a significant fraction of quasars (Population A; about 
50-60% of low z quasars). Our VLT spectra confirm that 
any dep endence of FWHM(H /3Br) on source luminosity 
is weak IjSulentic et inmOOi). Table El suggests that low 
redshift trends for asymmetries and line shifts are pre- 
served in the intermediate z sample. It is therefore not 
certain that FWIIM(H/3bc) is a valid estimator of the 
virial velocity for all sources even after proper H/?nCj 
[Oiii]AA4959,5007, Feii subtraction. 

A more physical approach to ATbh estimation uses 
the FWHM of the v ariable part of the H/3bc profile 
llPeterson et all l2004l) Fig , (l plots FWEM^ms values 
from IPeterson et alJ ^|2004^ versus our FWHM(H^bc) 
measures ( Marzia ni et al.L l2003a'l for all sources in com- 
mon. Sources with FWHM ^ 4000 km s~^ show a 
correlation while the situation is less clear for sources 
with FWHM^ 4000 km s"^ We recafl that this 
value indicates the nomin al population A-B bound- 
ary (jSulentic et al 1 l2000hl) that emerged in our El 
studies. A least-square best fit analysis yields a "cor- 
rected" FWHM(H/3bc) estimate: FWHM™"(H/3bc)~ 
-710(±800) + 1.13(±0.28) FWHM(H/3bc) for FWHM 
< 4000 km s-i and FWHM'=°"(H/3bc)~ 650(±1000) + 
0.79(±0.14) FWHM(H^bc) for FWHM ^ 4000 km s'^ . 
The slope depends somewhat on the fitting method (a ro- 
bust fit yields 0.94 and 0.72) but is always less steep for 
FWHM ^ 4000 km s~^. An immediate implication is that 
optically thick BLR gas responding to continuum changes 
shows a velocity dispersion correlated with - but slightly 
lower than that of the integrated profile at all FWHM. 

A break in the linear fit at 4000 km s~^ is consis- 
tent with several previous findings: (a) mean and pos- 
sibly systematic H/3rc profile differenc es between Pop. 
A and B sources (jSulentic et all l2002|) : ( b) the lack of 
stron g profile changes in Pop. A sources l)Giveon et all 



BH 



5.48-10^ 



ALa (5100A) 
10*4 ergs- 



H 0.67 



/FWHM(H^Bc) 
V 1000km s-i 



(7) 



where A/bh is in solar mass units (M0). 

Fig.Elshows the distribution of A/bh as a function of z. 
A/bh estimates are based on the FWHM(H/3bc) measures 



B sources. They might represent a distinct redshifted emis- 
sion component (Sulentic et al.. J002Jwhich may arise m 
less optically thick gas th an the rest of the H/3bc pro- 
file feule ntic et allEoO OcV If the redshifted component is 
emitted in an innermost VBLR, the FWHM of the whole 
H/3bc profile is obviously increased over the value due to 
the line component that is actually responding, and that 
is most likely located farther away from the central con- 
tinuum source. Even if the virial assumption holds for the 
VBLR (but the frequent asymmetries warn us that this 
might not be the case), the use of the FWHM from the 
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whole H/3bc profile and of tblr from the reverberating 
part, yields an Mbh overestimate. This interpretation of 
the correlations in Fig.^is advanced with caution because 
of the small sample size and especially because of the poor 
statistics for sources with FWHM > 4000 km s^^. It is 
also possible to produce a meaningful fit with a second- 
order polynomial. We apply a tentative correction to the 
FWHM measures FWHM™" (H/?bc) and therefore to re- 
sultant A/bh estimates using the second-order fit shown 
in Fig. 0] which approximates very well the linear trends. 

If an optically thin/no nvirialized component is present 
in the H/3bc profile of many sources then we might 
use as a virial estimator the FWHM of a line, or lines, 
arising in BLR gas but not likely to be present in the 
VBLR region. FWHM(FeiiA4570) is an obvious alter- 
native because there is no evidence for a VBLR emis- 
sion component in the broad Fell blends. It is not 
strictly correct to use it because the tblr- L\ rela- 
tion was deduced for H/3. The most serious difficulty 
lies in obtaining a rehable FWHM(FeiiA4570) estimate 
from the heavily blended Fell emission. Considering the 
EW and FWHM limits for detection of Fe iiopt emis- 
sion ijMarziani et al.l l2003a(l we conclude that a reason- 
able FWHM measurement is possible for «120 sources 
in our low z spectral atlas. Measurement uncertainties 
for FWHM(Feii) will be larger than for FWHM(H/5bc) 
and are estimated to lie between 20-50%. The best 
fit of FWHM(FeiiA4570) vs. FWHM(H^bc) (km s'l) 
is consistent with FWHM(FellA4570) « FWHM(H/5bc) 
for FWHM(H^Bc) - 4000 km s^^, while it is 
FWHM(FeiiA4570) « 0.67 • FWHM(H/3bc) + 820 km s'^ 
if FWHM(H/3bc) - 4000 km s^^. 

The relationship between FWHM(H^bc) and FWHM 
(Fell) confirms that individual Fell lines show approxi- 
mately the same width as H/3bc and as the rms H/3bc 
component implying a common kinematic environment if 
FWHM(H/3bc) ^ 4000 km s'^ . If FWHM(H/3bc) ^ 4000 
km s^^, FWHM(Feii) follows a trend closer to that of the 
rms H/3bc component. Therefore if the rms H/3bc compo- 
nent arises from gas in virialized motion, the same can be 
reasonably assumed for the whole FeiiA4570 emission. 

5.3. Mbh 3nd L/L-^^^ Dependence on Redshift 
5.3.1. Mbh 

Fig. 131 shows the distribution with z of Mbh esti- 
mates derived using uncorrected FWHM(H/?bc) mea- 
sures. This can be compared with Fig. where we 
show corresponding distributions of Mbh derived us- 
ing corrected FWHM'^°''''(H/3bc) measures and using 
FWHM(FeiiA4570). FWHM'=°"(H/3bc) values were de- 
rived from the second-order relation in Fig. ^ 

The distribution of data points in the log Mbh vs. z 
plane refiects selection efi'ects intrinsic to any fiux limited 
sample. If we consider AGNs that are radiating at a given 
L/LEdd we easily obtain Mbh for a given apparent mag- 
nitude: 



Mbh,90 = 3.85 • lO^h-^ (t^) ■ lO-"-*'"^ • 

\ LEdd J 

■ [1.5 {l~e-^) + (l-e"T3B6)]^ . + 

Appendix details how Eqn. |S1 was derived. If we 
adopt a limiting magnitude ttib ~ 17.5 (appropriate 
for the Hamburg-ESO survey) and consider a maximum 
Eddington ratio L/LEdd~ 1 we obtain the minimum Mbh 
detectable as a function of z. Limiting Mbk curves are 
shown in Fig. 13 Fig. |5l and Fig. El 

All sources show Mbh ^ 10^ Mq. In the range 1 ^ 
z ^ 2.5 we begin to find a significant black hole popu- 
lation in the range - lO^-^ M© < Mbh ^ 10^° Mq. If 
we consider Mbh values derived from FeiiA4570 as well as 
from FWHM™"(H/3bc) (Fig. E open and filled symbols 
respectively) we get a somewhat different picture. Almost 
all sources lie below « 5 • 10^ Mq (all sources with z ^ 1); 
only three sources whose Mbh has been computed with 
from FWHM'^°'''^(H/3bc) lie significantly above this hmit. 
Our sample of intermediate z quasars with compatible 
quality data is still small so caution is needed. It should 
be pointed out that the use of UPbc"^"" is self-consistent, 
since no H/3 surrogate line was used. 

The upper panel of Fig. identifies sources on the ba- 
sis of radio-loudness where open and solid symbols de note 
radio-loud (RL; see definition in lSulentic et allboO.*^ and 
radio-quiet (RQ) sources respectively. At low z RQ sources 
are distributed across the full range of Mbh up to Mbh^ 
10^ M(n, a lthough RL sources tend to have larg er Mbh 
ijMetcalf fc Magliocchetti . 2006; Marziani et all ,2003b, 
and references therein). The same consideration apply to 
Pop. A and Pop. B sources (lower panel of Fig. EJ, and 
this is not surprising since most Pop. A sources are RQ. 
At intermediate z, selection effects limit the detectable 
black holes to Mbh^ 10^ M©, so that only the high end 
of the quasar Mbh distribution can be traced. It is too 
early to decide whether systematic differences in the Mbh 
distribution of RQ and RL sources may still exist at in- 
termediate z, at least from the present data. 

5.3.2. L/LEdd 

Fig. 13 shows the distribution of L/Lsdd estimates as a 
function of z. A preliminary analysis in terms of L/LEdd) 
and including discussi on of bolometric lumino sity estima- 
tion, was presented in lMarziani et all l)2003b|) . 

The existence of apparently super-Eddington radia- 
tors in the low-z part of Fig. [3 involves the most extreme 
NLSyl sources. One caveat about interpreting any of them 
as super-Eddington involves the significant uncertainties 
associated with these estimates. If we assume a typical un- 
certainty of ±50% for the bolometric correction (neglect- 
ing beaming or lensing) and ±10% for the virial estimator 
(basically the uncertainty in the FWHM measurement) we 
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find a AlogL/iyEdd~ 0.13 (neglecting the scatter in Eqn. 

A serious source of uncertainty ignored until now in- 
volves source orientation with respect to the line-of-sight. 
The virial velocity dispersion that H/3bc is assumed to 
measure is now widely assumed to involve Keplerian ro- 
tation in an accretion disk or in a flattened distribution 
near the disk. If our candidate super-Eddington sources 
are mostly viewed with (face-on disk and/or pole-on jet) 
then Mbh can be significantly underestimated. An orien- 
tation correction applied to extreme NLSyl sources (we 
call them Pop. A "blue outliers" ) interpreted as face-on , 
will move them below L/L^m « 1 ijMarziani et al.lf2003bl 
see arrows in their Figure 12). At high z we so far find 
no strong evidence for any super-Eddington sources. The 
thick solid line in Fig. [T] shows the expected L/iEdd de- 
tection limit deduced from Eqn.|Slfor a flux- limited survey 
(ms w 17.5) and a quasar with Mbh w 4 • 10^ Mq. AU 
quasars in our sample fall above the minimum detectable 
L/L-Edd, suggesting that masses larger than « 4 • lO^M© 
are not strictly necessary from our data (see t ifi.ll for fur- 
ther discussion). 

A physic al basis for t he dichotom y between Pop. A 
and P op. B l)Sulentic et a l.. 2000b; M arziani et all l200ll 
[2003b) is supported by this analysis in the sense that Fig. 
^ shows evidence for a change at about FWHM(H/3bc) ~ 
4000 km s^^. Fig. 13 identifies Pop. A and B sources as 
filled and open symbols respectively. Here we apply the 
luminosity-dependent defi nition of the Pop. A -B bound- 
ary as defined in Fig. 6 of ISulentic et all l)2004|) . It shows 
that the Eddington ratio of Pop. A sources is systemat- 
ically larger than that of Pop. B, and that the apparent 
boundary between the two populations may increase with 
redshift. If we focus on the intermediate z quasars then 
Pop. A sources show an average L/Lsdd ~ 0.78 compared 
to 0.27 for 10 Pop. B sources. Even these small samples of 
Pop. A and B sources show significantly different L/L^dd 
distributions according to a Kolmogorov-Smirnov (K-S) 
test. RL AGNs are systematically low L/L^dd radiators 
since they are almost entirely Pop. B sources. 

It is interesting to note that with: (1) a limiting magni- 
tude TOB « 16.5, (2) L/LEdd^ 1 and (3) A/bh ^ 4-109 Mq 
we should detect fewer sources beyond z ~ 2 (no source 
below the dot-dashed curve of Fig.[7|). If our assumptions 
about the quasar bolometric correction are valid up to 
that redshift, selection effects may influence the relative 
frequency of Pop. A and B sources (the two Populations 
have different L/L^dd distributions) rather than the in- 
trinsic properties of LILs. In this case selection effects on 
L/LEdd should strongly influence the so-c alled "Baldwin 
effect " involving CivA1549 and other HILs l|Bachev et all 
12004 ^. and references therein) because CivA1549 is more 
prominent in Pop. B, and Pop. B sources are more easily 
lost at high z. 

6. Discussion 

The present paper provides Mbh estimates that have three 
advantages: (1) a consistent data analysis procedure is em- 



ployed over the entire redshift range 0.0 ^ z ^ 2.5 by 
using the same Mbh tracer, H/3bc; (2) S/N and resolu- 
tion of the spectroscopic data sample are high enough to 
permit a careful study of the H/3bc profile, and (3) the 
data quality allows reasonable estimates of FWHM Fell 
in many of the sources. 

6.1. What are the largest black hole masses? 



iNetzeJ l|2003l) discussed several problems with Mbh ~ 
10^° Mm If the black hole rnass vs . bulge mass (Mbuige) 
relation ((Ferrarese fc MerrittL I200(t1) is valid at high z then 
A^BH ^ 10^° Mq would impl y stellar velocity dispe rsion 
700 km s-i f following Icebhardt et aU l2000|) and 



resultant bulge m asses MhnVs ^ 10_ _Mf:^ which are not 
observ ed at low- z ijMcLure fc JarvisL l2004l iNetzeri 120031: 
Wan3, HoO^. Recent results for the fundamental plane 
of elliptical galaxies, and the most massive spheroids at 
z ^ 0.3 from Sloan Digital Sk y Survey measures, confirm 
that < 500 km s'^ fSDSS. iBernardi et all liooi 120051 
all galaxies with a ^ 500 km are likely due to chance 
superposition). 

There are several proposed interpretations of this prob- 
lem: (1) the Mbh ~ .^buige relationship may not strictly 
hold for all hosts, (2) the virial assumption is not appli- 
cable, (3) results are plagued by such large uncertainties 
- including the one of the lumino sity index a - that very 
large mass estimates are not real ()VestergaardL l2002fl , and 
(4) some systematic effects may not have been considered. 



6.1.1. The Mbh - Mbuigo Relationship 

Nuclei with A/bh 5 • 10^ M0 are not observed in galax- 
ies of the local Universe if a direct black hole mass de- 
termination is po ssible from circum-nuclear kinematics 
((Marconi fc Hunt] . i2Q(22.). They are expected to be rare 
and difficult to find, considering also that they should be 
in a dormant or nearly dormant stage at the present epoch. 
Integrating the qu asar luminosity function at z « 1.5 
l(Bovle et al.[l20od) . we find a comoving density of quasars 
above the HE limiting magnitude (which corresponds to 
AfB« -27.1 for a /c-correction a = 0.6) ~ 3 • 10"^ 
Mpc^"^. This indicates that the present-day density of 
the most massive black holes that were once luminous 
quasars should be very low, ~ 2 • 10"^ Mpc""^. Even so, 
they should be much more frequent than the very mas- 
sive spheroids that would host them if the A/bh^ Afbuigc 
relationship is valid. To estimate the density of spheroids 
with a^, ^ 500 km s^^, one must extrapo l ate th e a distri- 
bution function provided by'Sheth et al. (20 03l) to unob- 
served domains in a. Integrating the Sheth et al.l l(2003f) 
function for cr* ^ 500 km s~^, one obtains that the co- 
moving density of all local spheroids is three orders of 
magnitude lower (^ 6 • 10~^^ Mpc^'^) than that of the 
most massive black holes. A possible implication is that 
the A/bh ^ Afbuigc relation is not linear or of universal 
validity, i.e. some galaxies host larger black hole masses 
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than expected. This conclusion should remain valid also 
at intermediate redshifts. A constant MeH/Afbuige ratio 
cannot hold forever if bulge mass grows by secular evolu- 
tionary processes. Evolution goes in the sense that the 
spheroids m ay increase their masses at a later cosmic 
age {z <, 2; iTreu et all l2004t IPeng et all l2005(l signifi- 
cantly more than their central black hole: the most mas- 
sive bla ck hole should be al ready "brightly shining" at 
z « 2 ((McLure et al.L l200,4 and references therei n). A 
nonli near Mbh - -Mbuige relation has been proposed l|Laoii 
l200l|) for local hosts, with Mbh accounting for only 0.05 
% of the bulge mass in low-luminosity galaxies and 0.5 % 
in giant eUipticals. In this case, Mbh~ 10^° M0 would 
imply Mbuigo'^^ 10^^ Mq, which is at the upper end of 
the s pheroid masses measured locally l|Marconi fc Hunj . 

iool. 

The problem can be made to disappear if Mbh^ 3 • 

500 



ine proDiem can oe maae to aisappear it ai-q 
10^ M^. lBernardi et all ll2005l) find a maximum cr^ 
km s~^. These local galaxies (with expected Mbh ~ 3-10^ 
Mq) may have been the hosts of the luminous quasars 



at 1 < 
cr^f ^ 350 km s~^ is 



2 since the number density of spheroids with 
4-10~^ Mpc~^, still somewhat larger 
than the present-day density of supermassive black holes 
that were radiating at Mb ^ -25 10"'^ Mpc"'^). 

The curves in Fig. where Mbh^ 4 • 10^ Mq has 
been assumed, suggest that there is no need for Mbh ^ 
10^° Mq at z < 2.5. Similarly the brightest HE sources m 



the redshift range 1 



2 (mB ~ 15.0 at 2: w 1.7) are 



consistent with _L/LEdd~ 1 if Mbh~ 4 • 10^ Mq. In other 
words, if AIbh^ 10^° Mq, then a source at L/LEdd~ 1 
would be brighter than the brightest quasars observed in 



the range 1 ^ 
in this paper). 



2 (if Ho is close to the value assumed 



6.1.3. Statistical Errors 

It is interesting to consider in more detail the sources 
that show Mbh 10^° Mq in our sample. HE 0248-3628 
with the largest estimated Mbh is also one of the highest 
L/LEdd radiators with L/LEdd~ 1- It shows an anoma- 
lously large (by a factor ~10) flux with respect to our other 
VLT sources w hich makes it an extre mely luminous quasar 
(Mb^ -30; see lSulentic et al.Ll2004l) . It falls near the Pop. 
A-B boundary in Eigenvector 1 sp ace and is a borderliii e 
radio-loud source by our definition | Sulentic et all l2003l) . 
It disappears as an outlier if the Mbh estimated from the 
tabulated Mb is adopted. HE 2355-4621 is a second source 
with Mbh~ 10^° Mq. It behaves like a normal radio-quiet 
Pop. B source. The H/3bc profile shows a prominent red- 
ward asymmetry, which is strongly affecting the width at 
half maximum. It is also interesting to consider HE 1104- 
1805, which is the most luminous source in our new sam- 
ple with (Mb~ -29.5). It behaves like an ordinary Pop. A 
source and shows Mbh~ 5 • 10^ Mq which may be over- 
estimated if the continuum is lens brightened (Appendix 

El). 

The high-mass wing of the Mbh distribution in the 
redshift range 1 ^ z 2 is consistent with the wing of a 
Gaussian peaked at ~ 3 TO^ Mq and dispersion A log 
Mbh « 0.3 if Mbh is computed from FWHM™"(H/3bc) 
(4 -10^ Mq if no correction to H/3bc is applied). The 
Gaussian dispersion is consistent with the estimated er- 
rors of individual measurements. K — S tests do not favor 
significantly different peak masses or a much different dis- 
persion. This suggests that sources with Mbh ^ 3 -10^ 
Mq of our data might be mostly due to random errors 
associated to the uncertainty in individual Mbh measure- 
ments. 



6.1.2. Could A Non-virial Component Yield a Huge 

Mbh? 

One might envision a non-virial component that increases 
with source luminosity and systematically broadens the 
H/3bc profile. In our so-called Pop. A sources, the high 
ionization wind that can dominate CivA1549 emission 
might reasonably be expected to produce Balmer line 
emission as well. Such an additional (blueshifted) com- 
ponen t on the H/3Rn profile would increase FWHM mea- 
sures ((Sulentic et all .2006. HE 1249-0648 and HE 1258- 
0823 may show this effect). One intriguing possibility is 
that high-luminosity sources are extreme Pop. A sources 
which is supported by models invoking a radiation pres- 
sure driven wind. Such an effect would make the correc- 
tion (~ 500 km s"^) deduced for low-z Pop. A sources 
inadequate. At the other extreme, strong redward asym- 
metries observed in Pop. B sources indicate that the in- 
tegrated profile may be affec ted by gravitatio nal redshift 
and non-virial motions (Marz iani et alJ . f2003b|) . A quanti- 
tative analysis of these suggestions needs careful analysis 
beyond the scope of this paper, but in both cases any cor- 
rection would lower Mbh- 



6.1.4. Systematic Effects 

A first systematic effect considers the uncertainty in the 
index a. Our A^bh estimates show the onset of Mbh~ 
3T0^ Mq at z at z ^0.8 we observe some of the most 

luminous HE quas ars. The high luminosity range of the 
iKaspi et al.l l|2005l) relationship remains poorly sampled: 
there are just 2-3 sources in the range 45 ^ logAL^ ^ 
46 (to which most of our intermediate z sources belong), 
creating a sample bias and making the correlation analysis 
intrinsically unstable. In addition, a recent reanalysis of 
the rBT.R- \L\ correlation suggests a value of a as low as 
0.5 l|Vestergaard fc Petersonll2006|) . If a is overestimated 
by Aa « 0.17 over a luminosity range of ~ 10, Mbh rnay 
be overestimated by A log Mbh~ 0.2 in our intermediate 
z sources. 

Pop. A sources show good evidence that LILs are emit- 
ted in a strongly flattened system, probably an accretion 
disk or gas co-planar with the disk. CivA1549 in these 
sources seems to be dominated by a wind component 
l|Bachev et all l2004l) . If these considerations apply also 
to Pop. B sources (but it is by no means clear, given the 
large FWHM(H/3bc) of sources believed to be observed 
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pole-on) and if the maximum angle between disk axis and 
line-of-sight is w 45°, a correction could imply a factor 
^ 2 systematic increase in Mbh- However, these consider- 
ation of orientation effects may not even reopen the prob- 
lem of very large Mbh values: if our suggestion of a maxi- 
mum Mbh~ 3-10^ M0 is appropriate, taking into account 
the systematic orientation effects would yield a maximum 
Mbh~ 6 • 10^ Mq which is still plausible. The result that 
L/LEdd^ 1 at z ^ 0.8 would be reinforced by systematic 
orientation effects. 



Summing up, our data suggest that very large masses 
~ 10^*^ Mq may be not be real, and may be predominantly 
due to statistical errors and emission line profile broaden- 
ing that is in part non-virial. The data presented in this 
paper are consistent with Mbh not exceeding 3-10^ Mq for 
our sources if the correction to FWHM(H/3bc) described 
earlier is applied. 



6.2. The Best Mbh Estimators 

LILs like MgiiA2800 and Fell may yield more reliable re- 
sults than H/3bc- Hi Balmer line emission can be substan- 
tial from gas in a variety of physical conditions. In Pop. B 
sources, a very broad component may increase the FWHM 
of the integrated H/9Rr profile (also mimicked by lo w S/N 
data; iMcIntosh etalTllQQgt IShemmer et~all l2004|) . This 
very broad component may be optically thin to the ion- 
izing continuum, and therefore non-responsive to contin- 
uum changes. The H/?bc profile of Pop. A sources may 
be affected by a high-ionization component mentioned 
earlier. Fell is thought to be emitted in a region very 
optically thic k to Lyman continuum which is probably 
photo ionized ijVestergaard fc Petersonll2005HWang et al.l 
as the part of H/3bc responding to continuum 
changes should be. It is not surprising that the rever- 
berating part of H/3bc £^nd Fell provides width estimates 
which are consistent, since they are expected to mea- 
sure the width of a similar sub-region with in the BLR . 
Simil ar considerations apply to MgilA2800 JWills et all 
Il985^ since Mgll A2800 should be mainly emitted in the 
same zone as Feii. lMcLure fc Dunlon (20041 present virial 
A/bh estimates for « 13000 quasars in the redshift inter- 
val 0.1 ^ z ^ 2.1 based on spectra from the SDSS first 
data release. The mean Mbh increases with increasing red- 
shift basically as shown in Fig. |21 The mass values found 
by them are also consistent with a limiting A/bh around 
3 • 10^ Mq, with large scatter. They use FWHM(H/3bc) 
or FWHM{MgllA2800) and find a consistency between the 
most massive at z « 2 and those at z « 0. Given measure- 
ment difficulties, and doubts about the virial assumption 
for most other lines, corrected measures for H/3bCi Fell, 
and MgiiA2800 may offer the best hope for reliable A/bh 
and L/L-Edd estimates out to z ^ 2.5. 



7. Conclusions 

Nine intermediate z VLT/ISAAC spectra with high res- 
olution and S/N supplement an earlier sample of 17 
sources. Emission line measurements on H/3bc and Fell 
presented in this p aper strengthen the conclusion of 
ISulentic et al that luminosity effects are weak or 

absent in the low- ionization lines of AGNs. Results on the 
H/3bc profile are consistent with the population A-B hy- 
pothesis and Eigenvector 1 parameter space concept de- 
veloped for low z AGNs. We computed virial masses and 
Eddington ratios for the 25 intermediate-z objects plus 
about 280 lower z sources, using the same emission line 
over the entire redshift range for the first time. 

We also have how the distributions of A/bh and 
L/L-Edd vs. z are shaped by selection effects intrinsic to 
any fiux-limited survey at the low Mbh end. At the high 
Mbh end, masses exceeding a few 10^ Mq may be rare 
if corrections for non-virial broadening and statistical er- 
rors are taken into account. This suggestion is based on 
just 25 objects distributed over the entire redshift range 
0.9 ^ z ^ 2.5. Confirmation from a larger sample of 
intermediate-z observations is needed. 
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ally lensed quasar with the Space Telescope Imaging 
Spectrograph (STIS). It is a doubly imaged QSO with 
a source redshift of z = 1.58 and an image separation 
of 0" .644. The flux ratio ^ of the lensed images shows a 
strong dependence on wavelength. In the R and I bands, 
A is brighter than B by about 0.45 mag while the two 
are almost equal in the band B. For smaller wavelengths, 
especially close to the limit near 2000 A, B becomes 
much brighter than A by about 1.3 mag. A natural ex- 
planation for this effect is differential reddening caused 
by different extinc tion effects in the two lines of sight 



pact ob ject in the lensing gala xy also plays a role in this 
source ijWucknitz et al.l 12003^ . Unfortunately, the small 
separation does not allow to us to distinguish component 
A and B on the ISAAC spectrum. The HE 0512-3329 ac- 
quisition image is compatible with an unresolved source. 
Considering the flux ratios in the R and I band, it is rea- 
sonable to assume that our spectrum is dominated by the 
A component. 

A. 2. HE1104-1805 

HE 110 4-1805 is a double- irn age lensed quasar discov- 



A9 = 3". 19, the source redshift is z,^ = 2.319 , and the 
lens redshift is zi = 0.729. IWisotzki et alJ l)l995 h reported 
that the continuum flux in both images is highly variable 
but that the line fluxes do not change, as expected if mi- 
crolensing is operating. On the acquisition image of HE 
1104 there is a second source at 3.5" but it is completely 





separation is 



off-slit. 
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(1) 
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(4) 


(5) 


(6) 
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(8) 


(9) 


do) 


fii) 
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HE0507- 


3236 


17.36 


1.5770 


(7) 


1 


—27.6 


^ 0.51 


2003-10-22 


J 


180 


18 


1.23-1.09 


25-20 


HE0512- 


-3329 


17.03 


1.5873 


(7) 


1 


-28.0 


<, 0.38 


2003-10-22 


J 


180 


16 


1.07-1.02 


20-15 


HE0926- 


-0201 




i.Doz4 




1 


on n 
— za.U 


^ -U.oo 




T 
J 


1 on 




111 117 
i.ii-i.i / 


QA 1 n 
oU— iU 


HE0946- 


-0500 


16.24 


1.1013 


(7) 


1 


-28.0 


<, 0.03 


2004-02-11 


sZ 


180 


12 


1.17-1.27 


30 


HE1003+0149 


16.45 


1.0809 


(7) 


1,2 


-27.7 


0.20 


2004-02-13 


sZ 


180 


12 


1.19-1.14 


20 


HE1017- 


-0009 


16.69 


1.1295 


(7) 


1 


-27.6 


0.87 


2004-02-28 


sZ 


180 


12 


1.97-1.62 


20-10 


HE1104- 


-1805 


16.45 


2.3192 


(7) 


1,2 


-29.6 


<, 0.19 


2004-02-11 


sH 


180 


20 


1.14-1.35 


20 


HE1249- 


-0648 


16.72 


1.1940 


(7) 


1 


-27.7 


^ 0.23 


2004-02-05 


sZ 


180 


16 


1.04-1.06 


30-15 


HE1258- 


-0823 


16.26 


1.1632 


(7) 


1 


-28.1 


<, 0.04 


2004-03-29 


sZ 


180 


17 


1.06-1.16 


30 



Apparent B magnitude corrected because of Galactic absorption. 
^ Redshift, with uncertainty in parenthesis. 
" Lines used for redshift calculations: 1: H/3, 2: [Olll]A5007. 

Absolute B magnitude, computed for _ffo=70 km s~^Mpc~^, Q.ia = 0.3, = 0.7, and fc-correction spectral index a=0.6. 

Decimal logarithm of the specific flux ratio at 6cm and 4400 A (effective wavelength of the B band). Upper limits are from 

the NVSS (« 2.5 mjy), and would place all undetected sources in the RQ domain. 
' Date refers to time at start of exposure. 
^ Photometric band. 

^ Detector Integration Time (DIT) of ISAAC, in seconds. 

' Number of exposures with single exposure time equal to DIT. 

' Airmass at start and end of exposure. 

S/N at continuum level in the proximity of H/3. Two values are reported in case of different S/N on the blue and red side 
of H/3 (blue side first). The S/N value is with N estimated at a 2(7 confidence level i.e., 2 times the rms. 



Table 2. Measurements of Fluxes, Equivalent Widths and FWHM of Strongest Lines 



Object name 


F(H/3bc)=' 




W(FenA4570)'^ 


FWHM(FenA4570)'' 


(1) 


(2) 


(3) 


(4) 


(5) 


HE0507-3236 


12.5 ± 3.0 


62 ±20 


17 ± 3 


4100 1 


HE0512-3329 


20.0 t'il 


7C -1-30 
'0 -10 


47 ± 10 


2400 


HE0926-0201 


26.0 ± 4.0 


72 ± 15 


20 ± 2 


3400 ± 1200 


HE0946-0500 


17.0 til 


59 tZ 


19 1^ 


2900 ± 1000 


HE1003+0149 


7.0 ± 1.0 


48 ± 10 


67 ± 12 


2100 ± 1400 


HE1017-0009 


10.5 ± 1.0 


60 


11 tl 


3600 t'^l 


HE1104-1805 


20.0 ± 4.0 


72 ± 15 


39 +1 


4100 ± 1700 


HE1249-0648 


15.5 ± 2.0 


63 ± 15 


31 tl 


1400 


HE1258-0823 


24.0 ± 3.0 


44 ± 10 


20 ± 4 


1400 tir 



Appendix B: Minimum Mbh a"d L/L^dd as a 
Function of z 

The absolute B magnitude Mb can be related to the spe- 
cific luminosity at 5100 A assuming an average spectral 
shape between 4400 (effective wavelength of B band) and 
5100 A. If the spectral s hape is described bv a p ower-law 
{U fx ly-'') with b = 0.3 ijMarziani et al.U2003b|) : 

log[ALA(A = 5100A)] = -0.4Mb + 35.497. (B.l) 

If we further assume that L ~ IDA • L\, with A=5100 
A i.e., a bolometric correction factor 10, we can write: 

log[ALA(A = 5100A)] = —4 1.3 • 1O47Mbh,90 (B.2) 

10 • L^dd 



where the black hole mass A/bh has been written in 
units of 10^ Mq. The absolute magnitude M-q is: 

Mb = me -t- 5 - SlogdL - 2.5(a - 1) log(l + z) (B.3) 

where the luminosity distance c^l is in parsec, and the 
last term is the K{z) correction. 

Assuming f^M 7^ 0, 17a 7^ 0, and flk = 0, we have for 
the co-moving distance; 

c f^' dz' , 
dc = -rr B.4 

This integral can be approximated with residuals less 
than 3% at all z (maximum error for z — > 0). 

do - ^ [(1.500) (1 - e--^) + 0.996 (l - e-i-He )1 (B.5) 
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Object name 


F(H/3nc)'= 


W(H/3nc)' 


W([Oiii]A4959)f 


W([Oiii]A5007)' 


(6) 


(7) 


(8) 


(9) 


(10) 


HE0507-3236 


7: 


0.3: 


2.7 ± 0.5 


8.4 ± 0.6 


HE0512-3329 


3: 


0.1: 


2.5 ± 0.5 


6.5 ± 1.0 


HE0926-0201 


28: 


0.8: 


5.6 ± 0.5 


15.2 ± 1.0 


HE0946-0500 


14 +1 


0.9 


1.1 ± 0.5 


3.4 ± 0.6 


HE1003+0149 


15 


1.0 


0.6 ± 0.5 


1.7 ± 0.5 


HE1017-0009 


15 +_f 


0.8 


1.4 ± 0.5 


3.8 ± 0.6 


HE1104-1805 


60 ± 18 


2.0 ± 0.5 


3.9 ± 0.5 


11.4 ± 0.5 


HE1249-0648 


18 ± 5 


0.7 ± 0.3 






HE1258-0823 


9: 


0.8 ± 0.3 







^ Rest frame flux of H/3bc in units of 10~^^ ergs s~^ cm~^ and ±2a confidence level uncertainty. 

^ Rest frame equivalent width of H/3bc in A ±2(t confidence level uncertainty. 

" Rest frame equivalent width of the FellA4570 blend in A ±2a confidence level uncertainty. 

FWHM of lines in the FellA4570 blend and uncertainty at 2a, in km s~^. See text for details. 
" Rest frame flux of H/3nc in units of 10~^^ ergs s~^ cm~^ and ±2ct confidence level uncertainty. Colons indicate highly 

uncertain values. 

' Rest frame equivalent width of H/3nc, [Oiii]14959, and [Olll]A5007 in A, with uncertainty at 2a. Colons indicate highly 
uncertain values. 



Table 3. H/3bc Line Profile Measurements 



Source 


FWZP 




FWHM=' 




A.I." 




Kurt.'* 




c(0/4) 


^a,b 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


HE0507- 


-3236 


16000 


2100 


3500 


300 


-0.15 


0.10 


0.31 


0.06 


-500 


1100 


HE0512- 


-3329 


16000 


3700 


3100 


300 


-0.16 


0.07 


0.31 


0.04 


-100 


1900 


HE0926- 


-0201 


17000 


1400 


5100 


500 


0.21 


+0.07 

-0.19 


0.27 


0.06 


1200 


700 


HE0946- 


-0500 


20000 


2800 


3600 


250 


-0.06 


0.09 


0.37 


0.05 


-600 


1400 


HE10034 


-0149 


9000 


2100 


2900 


280 


-0.06 


+0.07 
-0.11 


0.33 


0.05 


-400 


1000 


HE1017- 


-0009 


16000 


3800 


6200 


290 


-0.04 


+0.08 
-0.05 


0.41 


0.04 


200 


1900 


HE1104- 


-1805 


15000 


1500 


4300 


260 


0.11 


0.09 


0.38 


0.05 


900 


800 


HE1249- 


-0648 


21000 


1000 


4900 


500 


0.13 


+0.11 
-0.07 


0.31 


0.05 


1400 


500 


HE1258- 


-0823 


17000 


2600 


4400 


500 


-0.08 


+0.20 
-0.08 


0.23 


0.06 


-700 


1300 



The luminosity distance is: 
dL = dc(l + ^) (B.6) 

and we can use the new relationship to derive a working 
relationship for Mbh^ 

Mbh,90 = 3.85 • 10^/1-2 ("t^) ■ lO""-^"^ • 

• [1.5 (1 - e^^) + (1 - e'^)f (I -f z)^^- 

where h = Ho/75. Curves shown in Figs. El El El 13 
assume a — 0.6. Differences for 0.3 ^ a ^ 0.6 are minor. 
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-300 
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Fig. 1. Calibrated VLT-ISAAC spectra for 9 new intermediate- redshift quasars. Abscissae are rest-frame wavelengtli 
in A, ordinates the rest-frame specific flux in units of 10~^^ ergs s~^ cm~^ A~^. 
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Fig. 2. Spectral Atlas for the new intermediate-redshift quasars. Left-hand panels show the continuum-subtracted 
H/3 spectral region. Axes are the same as Fig. ^ The best-fit Fell emission model is traced as a thin (green) line. 
Right-hand panels show the enlarged H/3 profiles after continuum and Fell subtraction. The (blue and red) thick line 
shows a spline fit of the short and long wavelength sides of the H/?bc profile, respectively. 
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Fig. 3. Distribution of Mbh estimates as a function of z using uncorrected FWHM(H/3bc) values. Intermediate z 
sources observed with ISAAC are at z ^ 0.8. The soUd curve estimates the minimum detectable Mbh for our sample 
- see text. 
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Fig. 4. FWHM meas ured from the variable part of the H/3bc profile for low-z sources wit h reverberation data 
I Peterson et all EoO^ versus FWHM(H/3bc) (km s~^) measures from iMarziani et ^ l|2003cf) . Solid straight hnes 
show best fits for FWHM(H/3bc)< 4000 km s~^(Pop. A) and for sources with broader hnes(Pop. B). The solid curve 
shows a weighted least-square fit of all data with a second order polynomial. The diagonal dot-dashed line indicates 
the location of equal ordinate and abscissa values. 
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Fig. 5. A/bh computed from corrected FWHM'^°'''^(H/3bc) as a function of z (filled circles). A^bh estimates computed 
from FWHM(FeiiA4570) measures are shown as open symbols. Symbol types and solid curve are as described for Fig. 
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Fig. 6. A/bh computations using same samples as Fig. O where RQ and RL sources are indicated by filled and open 
symbols respectively (upper panel). Triangles indicate sources for which radio data are insufficient to compute a 
meaningful i?K- The lower panel distinguishes Mbh for Pop. A (filled symbols) and Pop. B (open symbols). The solid 
curve is as described for Fig. 13 




Fig. 7. Distribution of L/LEdd with z following Fig. |21 with XLx values derived from the 5100 A rest-frame flux and 
with a = 0.67. El Populations A and B are indicated by filled and open symbols respectively. The solid line traces 
the minimum L/L^dd for a fixed mass value of 4 • 10^ Mq to observe a quasar above the limiting magnitude of the HE 
quasar survey tob ~ 17.5. The dot-dashed lines indicate the same limit but assuming limiting magnitude tub w 16.5. 



